Virus Evolution

The Capripoxviruses (CaPV) comprise three species: goatpox virus (GTPV), sheeppox virus (SPPV) and lumpy skin disease virus (LSDV). They are large double-stranded DNA viruses with highly conserved core genomes and variable terminal regions. Previous studies have described variation in CaPV gene content, their broader population structure and the contribution of non-coding and structural variation remains opaque. This study investigated the genomic diversity and evolutionary history of GTPV and SPPV using an integrative framework combining phylogenetics, pangenome variation graphs (PVGs), and gene-specific analyses. We found marked differences in population structure between the two viruses. GTPV comprised three deeply divergent and genetically stable lineages with limited evidence of recent gene flow, whereas SPPV had weaker clade separation consistent with an ancestral bottleneck followed by recent population expansion. PVG-based analyses indicated that GTPV has a comparatively closed pangenome, while SPPV remains open, particularly at the genome termini. Structural and haplotypic variation was concentrated at the inverted terminal repeats (ITRs), which moderate host immunity and specificity. In several lineages, extended putative ORFs spanning adjacent terminal genes were observed, indicating recurrent structural plasticity at the genome ends. Patterns of gene-specific conservation and divergence highlighted loci under strong constraint and lineage-specific structural changes that may contribute to host specificity. Together, these results demonstrate how graph-based genome models complement gene-based analyses in resolving poxvirus genome evolution and provide a resource for improved comparative and population genomic studies of large DNA viruses. SignificanceCapripoxviruses are economically important livestock pathogens, yet the genomic mechanisms underlying their diversification and host specificity remain poorly resolved. By applying pangenome variation graphs alongside phylogenetic and gene-level analyses, this study reveals fundamental differences in how goatpox and sheeppox viruses have evolved. Goatpox virus had a deeper, more stable lineage structure, whereas sheeppox virus was more recent and diverse. Importantly, structural variation at the inverted terminal repeats emerged as a major driver of genomic diversity, including lineage-specific haplotypes and variable gene structures. These findings demonstrated the value of graph-based genome representations for resolving complex variation in large DNA viruses and provides a framework for improving genomic surveillance, comparative analyses, and future investigations into host range, virulence and tropism.

Understanding how viral communities vary across co-occurring hosts and environments is essential for assessing species-specific viral risks under changing land use and climate. This is particularly relevant for managing introduced bees, which face persistent viral threats themselves, as well as transmitting plant viruses. Here, we compare RNA viromes of the long-established honeybee (Apis mellifera, introduced to Tasmania in 1831) and the more recent invader, the bumblebee (Bombus terrestris, invasive since 1992), across 14 Tasmanian sites - an island still free of the viral vector, Varroa destructor. Using a metatranscriptomic approach on total RNA from whole bees, we identified insect- and plant-associated viruses and inferred phylogenetic patterns of insect viral sharing, divergence, and potential cross-species transmission. We also assessed spatial and environmental drivers of viral composition, diversity, and richness. Geographic longitude, precipitation, temperature, and pasture percentage influenced the total, insect-, and plant-associated viromes of B. terrestris. In contrast, for A. mellifera, only precipitation and temperature were associated with insect and plant viral alpha diversity and community composition. Phylogenetic analyses revealed that Black Queen Cell virus in A. mellifera from Tasmania has diverged from mainland Australian sequences, and two distinct sub-strains of Lake Sinai virus 1 were shared by both bee species. Lake Sinai virus 3 showed evidence of interspecies transmission between A. mellifera and B. terrestris. Notably, this study provides the first detection of Moku virus in Australian bees and globally in bumblebees, suggesting potential interspecies transmission among social Hymenoptera. Overall, our findings demonstrate local viral diversification and reveal that B. terrestris viromes are more strongly shaped by environmental factors than those of A. mellifera, underscoring the importance of monitoring invasive pollinators as reservoirs and vectors of viral emergence.

The antigenic evolution of human seasonal influenza viruses is primarily driven by single amino acid substitutions immediately adjacent to the receptor binding site in the hemagglutinin (HA) protein. The ability to predict these substitutions would allow vaccine strains to be selected with an understanding of likely future antigenic variation. Here, we estimate the effect of HA substitutions on viral fitness using measurements of convergent evolution in a large phylogeny. We show that the substitutions which have historically caused major antigenic changes in H3N2 influenza viruses were nearly always one of few substitutions near the HA receptor binding site estimated to be under positive selection in sequences collected before the antigenic transition, based on convergent acquisition of the substitution in multiple independent lineages. Furthermore, this signal predates the establishment of the major clade containing the antigenic substitution by more than one year, so is highly informative for prospective prediction.

In most paramyxoviruses, RNA editing in the P gene enables expression of the V protein. Human parainfluenza virus type 1 (HPIV-1) differs from most paramyxoviruses in that it lacks RNA editing and does not produce a functional V protein, although its genome retains sequences corresponding to the ancestral V reading frame. Here, we analyzed all HPIV-1 genome sequences available in the NCBI GenBank database to assess the evolutionary state of this V protein-specific region. Using Sendai virus (SeV) as a closely related reference with an identical P gene length, we defined a pseudo-V reading frame by virtually inserting a single nucleotide at the conserved RNA editing site. In this pseudo-V frame, HPIV-1 showed a marked excess of stop codons within the 253-amino-acid region corresponding to the post-editing sequence, far exceeding expectations under random codon usage. This pattern was not observed in other viral genes analyzed under the same definition, nor in SeV, nor was it reproduced by in silico evolutionary simulations under constraints preserving the primary open reading frame. These results are consistent with a virus-specific evolutionary trajectory following the loss of RNA editing, rather than with generic coding constraints acting on overlapping reading frames.

Understanding the dynamics of HIV epidemics is important to control them effectively. Classical methods that mainly rely on occurrence data are limited by the fact that an unknown part of the epidemic eludes sampling. Since the early 2000s, phylodynamic methods have enabled the estimation of key epidemiological parameters from virus genetic sequence data. These methods have the advantage of being less sensitive to partial sampling and to provide insights about epidemic history that even predates the first samples. In this study, we analysed 2,205 HIV sequences from the French ANRS PRIMO C06 cohort. We identified and were able to reconstruct the temporal dynamics of two large clades that represent the HIV-1 epidemics in the country. Using Bayesian phylodynamic inference models, we found that the first clade, from subtype B, originated in the end of 1970s, grew rapidly during the 80s before decreasing from 2000 to 2015 and stagnating since then. The second clade, from circulating recombinant form CRF02_AG, emerged and spread in the 80s, grew again in the early 2000s, before declining slightly. We also estimated key epidemiological parameters associated with each clade. Finally, using numerical simulations, we investigated prospective scenarios and assessed the possibility to meet the 2030 UNAIDS targets. This is one of the rare studies to analyse the HIV epidemic in France using molecular epidemiology methods. It highlights the value of routine HIV sequence data for studying past epidemic trends or designing public health policies.

We used 2492 whole genome sequences from Laos to investigate the molecular epidemiology of SARS-CoV-2 from 2021 through 2024, covering the major waves of COVID-19 disease in Laos including time periods of travel restrictions and after relaxation of travel across international borders. We identify successive waves of COVID-19 caused by shifts in the dominant lineage, beginning with the Alpha variant in April 2021 and continuing through the Delta and Omicron variants. We quantify a shift from a small number of viral introductions responsible for widespread transmission in early waves to a larger number of introductions for each variant after travel restrictions were lifted, and identify potential routes of introduction into the country. Our study underscores the importance of genomic surveillance to public health responses to characterize viral transmission dynamics during pandemics.

Prior to 2024, highly pathogenic avian influenza H5N1 clade 2.3.4.4b viruses circulated predominantly in wild birds and poultry. In 2024 and 2025, 2.3.4.4b genotypes B3.13 and D1.1 were detected in United States dairy cattle. Using whole-genome and segment-specific phylodynamic inference, we estimate that B3.13 and D1.1 spilled over from wild birds into dairy cattle in late 2023 and late 2024, respectively. Spillover occurred shortly after the formation of the reassortant genotypes and was followed by months of cryptic transmission prior to detection. We found that both B3.13 and D1.1 evolved at higher rates in cattle relative to birds, primarily due to relaxed purifying selection. Site-specific analyses identified genomic sites under positive selection in cattle relative to birds, indicating adaptation and likely contributing to improved viral fitness after spillover. Intensified genomic surveillance in dairy cattle is essential as population immunity introduces additional selection pressures, with ever-changing risk for human emergence.

Sliding-window phylogenetic analyses of multiple sequence alignments (MSAs) generate sequences of phylogenetic trees that can reveal recombination and other sources of phylogenetic conflict, yet comparing trees across genomic windows remains challenging. Phylo-Movies is a browser-based tool, also available as a standalone desktop application, that decomposes topological differences between consecutive phylogenetic trees into interpretable subtree migrations and animates these transformations. We demonstrate its utility in two contexts: identifying recombintion breakpoints in norovirus genomes, where lineages shift from polymerase-based to capsid-based clustering at the ORF1/ORF2 junction, and detecting rogue taxa that change position across bootstrap replicates. Phylo-Movies complements summary statistics such as Robinson-Foulds distances by showing which lineages move, where they move from, and which new groupings they form. Phylo-Movies is freely available at https://github.com/enesberksakalli/phylo-movies, with a norovirus demonstration video at https://vimeo.com/1162400544, the first rogue taxon example at https://vimeo.com/1162561152, and the second example at https://vimeo.com/1162563101.

The yerba mate psyllid (Gyropsylla spegazziniana) poses a significant threat to yerba mate crops, causing extensive economic losses. While some ecological aspects as well as control strategies have been studied, its associated viral diversity remains unexplored. Here, by generating the first RNA high-throughput analysis (HTS) of this pest, we explored the G. spegazziniana virome, revealing novel and diverse RNA viruses. We characterized five new viral members belonging to distinct families, with evolutionary cues of beny-like viruses (Benyviridae), picorna-like viruses (Picornaviridae), and sobemo-like viruses (Solemoviridae); which were tentatively named Gyropsylla spegazziniana beny-like virus 1 (GSBlV1), Gyropsylla spegazziniana picorna-like virus 1 (GSPlV1), and Gyropsylla spegazziniana sobemo-like virus 1-3 (GSSlV1-3), respectively. Phylogenetic analysis of the bi-segmented and highly divergent sobemo-like viruses showed a distinctive evolutionary trajectory of its encoding proteins at the periphery of recently reported invertebrate Sobelivirales. The beny-like virus belonged to a cluster of insect-associated beny-like viruse; while the picorna-like virus clustered together with psyllid-associated picorna-like viruses. These results highlight the existence of a complex virome within G. spegazziniana and establish the basis for future studies investigating the ecological roles, evolutionary dynamics, and potential biocontrol applications of these viruses in the G. spegazziniana -yerba mate eco-systems.

Sarbecoviruses, including SARS-CoV and SARS-CoV-2, are frequently linked to Rhinolophus bats as their putative natural reservoirs. Angiotensin-converting enzyme 2 (ACE2), a host carboxypeptidase widely expressed in mammalian tissues, plays a critical role in sarbecovirus infection by serving as the cellular receptor for the viral spike (S) protein. Given recent human outbreaks and pandemics caused by members of sarbecoviruses, and the wide distribution of Rhinolophus bats, it is essential to maintain surveillance of these viruses while improving our understanding of their interactions with bat hosts, particularly the ACE2 receptor. However, while Rhinolophus bats from Asia have been relatively well studied, African Rhinolophus bats remain underrepresented and require further investigation. In this study, five Rhinolophus bat lung samples were obtained from Zambia, and ACE2 genes from these individuals were cloned and sequenced. We further evaluated the susceptibility of ACE2 variants to a panel of sarbecoviruses, revealing key residues that influence viral infectivity. ACE2 polymorphism was observed among Rhinolophus simulator individuals, revealing multiple ACE2 genotypes within the sampled population. However, R. simulator ACE2s did not permit infection by the clade 3 Afro-Eurasian sarbecoviruses tested in this study. Notably, RhGB01 and BM48-31 virus utilized only Rhinolophus blasii ACE2. Mutational analyses further suggested that ACE2 residues 31 and 41 play important roles in modulating spike-ACE2 interactions. This study reports 4 unique ACE2 sequences of R. simulator and R. blasii, and provides new insights into the molecular interactions between African Rhinolophus species ACE2s and the S protein of sarbecoviruses circulating in Africa and Europe. ImportanceAs putative natural reservoirs of sarbecoviruses, including SARS-CoV and SARS-CoV-2, Rhinolophus bats play a critical role in the emergence of zoonotic coronaviruses, making it essential to understand their interactions with these viruses for future pandemic preparedness. While Asian Rhinolophus bats have been relatively well studied, African species remain underrepresented, highlighting the need for further investigation. In this study, we cloned and sequenced ACE2 genes of five Rhinolophus bats collected in Zambia, Africa. We identified ACE2 polymorphism among Rhinolophus simulator individuals, although this variation was not associated with susceptibility to the clade 3 Afro-Eurasian sarbecoviruses examined. In addition, we identified key ACE2 residues that govern SARS-CoV-2 spike-ACE2 interactions and contribute to distinct infectivity patterns across species. These findings expand our understanding of the molecular determinants of sarbecovirus host range and support improved surveillance and risk assessment of emerging coronaviruses.

Oropouche virus (OROV), a member of the Peribunyaviridae family endemic to South America, is a current public health threat. The recent OROV outbreak driven by a novel reassortant has caused a dramatic increase in cases in 2024 (13,785 in Brazil, versus only 261 from 2015-2022) with sustained levels of transmission in 2025. Previously underreported outcomes have been recognized including miscarriage, microcephaly, encephalitis, and death. OROV lethality in humans has been attributed to severe coagulopathy with liver involvement, and epidemiological data suggests acute hepatitis occurs in mild cases of Oropouche fever, highlighting the underrecognized role of the liver in OROV pathogenesis. We present two discrete mouse models of OROV hepatic disease -- a lethal model that recapitulates the severe coagulopathy seen in fatal human cases and a model of self-resolving acute hepatitis which recapitulates mild human disease. OROV causes focal hepatic necrosis in mice, which progresses to massive necrosis and death when the Type I interferon receptor is antagonized. Additionally, we found a contemporary OROV isolate is less pathogenic in mice than a historic prototypical strain. These studies enhance our understanding of OROV pathogenesis and pave the way for potential therapeutic development and evaluation. ImportanceThe disease burden of Oropouche fever has been underrecognized and underreported, as highlighted by the increased testing seen in the ongoing outbreak. Specifically, the role of the liver in Oropouche virus pathogenesis has been neglected. Given the ongoing outbreak and increase in severe disease manifestations, there is a present need to understand Oropouche virus pathogenesis and test potential therapeutics. The mouse models of Oropouche-induced hepatitis presented here provide a means to understand how Oropouche virus causes liver damage in both a lethal and sublethal context. These models will be useful for the preclinical evaluation of vaccines and therapeutic treatments. Additionally, we compare the pathogenicity of a historical Oropouche virus isolate to a contemporary human isolate in a lethal mouse model. This represents an additional step towards understanding whether the circulating Oropouche virus isolates are uniquely more pathogenic or if increased testing has highlighted previously unreported outcomes.

The COVID-19 pandemic saw successive emergence and global spread of novel viral variants, exhibiting enhanced transmissibility or evasion of immunity. While the genotypic and phenotypic basis of SARS-CoV-2 variants have been extensively characterized, the evolutionary factors governing their patterns of emergence are less well understood. In this study we systematically investigated how the invasion dynamics of viral variants depend on variant phenotype (increased transmissibility or immune evasion), source (local evolution vs importation), the timing of introduction, the distribution of population susceptibility, and the contact network structure. Using a stochastic multi-strain epidemic model, we find that strains with only a transmission advantage are more likely to emerge earlier in the epidemic, and rapidly and predictably dominate the viral population. In contrast, immune-escape variants tend to linger at low prevalence for extended time periods after emergence, avoiding detection, until a critical amount of immunity has built up in the population and they begin to rapidly outcompete existing strains. We find that two common features of realistic human contact networks---heterogeneity in contacts (overdispersion) and clustering---lead to more punctuated evolutionary dynamics. This work provides insight into past dynamics of SARS-CoV-2 variants and can help define planning scenarios for future epidemic modeling efforts.

Understanding how genetic variation in pathogens influences clinical phenotypes observed in infected hosts is a fundamental challenge in evolutionary genomics and public health. Phenotypic traits such as infection severity are often non-randomly distributed within the pathogens phylogeny, suggesting the existence of evolutionary determinants but also violating the independence assumption underlying classical genome-wide association studies and potentially leading to inflated false positive rates. We present MutaPhy, a phylogeny-based method aimed at detecting correlations between a binary host phenotype and the corresponding pathogen genome by directly utilizing the hierarchical structure of phylogenetic trees. MutaPhy encompasses three different scales: (i) a subtree scale, on which relevant clades over-representing the phenotype of interest are detected using permutation-based tests; (ii) a tree scale, which agglomerates local signals into a global association statistics; and (iii) a site scale, whereby candidate mutational events on branches leading to significant clades are examined using ancestral sequence reconstruction. We evaluate the statistical behavior and detection performance of MutaPhy using simulations under diverse evolutionary scenarios. We also compare this tool to several existing phylogenetic association methods. As illustrative applications, we apply MutaPhy to dengue virus and hepatitis C virus datasets associated to clinical phenotypes in human hosts. Our results highlight the ability of the proposed approach to detect viral lineages associated to over-represented phenotypes while revealing limited evidence for robust mutation-level associations in these particular datasets. Altogether, MutaPhy provides a framework for guiding genotype-phenotype association analyses by leveraging phylogenetic structure, thereby reducing false positive findings and improving the interpretability of association signals.

The emergence of the panzootic clade of highly pathogenic avian influenza H5N1 (2.3.4.4b) in 2020 marked a major expansion in the host range of influenza A viruses (IAVs), raising concerns about further cross-species transmission events and zoonotic spillover. Introduction of 2.3.4.4b viruses into U.S. dairy herds has resulted in widespread circulation, accompanied by reduced milk yield, mastitis, and high viral loads in milk. Notably, virus circulation in dairy cattle represents a novel route for mammalian adaptation and transmission that has already led to more than 40 human cases in the U.S. since 2024. Here, we investigated whether avian clade 2.3.4.4b viruses could infect mammary tissue from Aberdeen Angus, Holstein Friesian, and Limousin cattle, three breeds commonly farmed in Europe, the Americas, and Oceania. Using mammary gland explants, we inoculated tissues with attenuated reassortant viruses expressing the haemagglutinin and neuraminidase glycoproteins of three 2.3.4.4b viruses that predated the emergence of H5N1 in US cattle: A/chicken/England/053052/2021 (AIV07), A/chicken/Scotland/054477/2021 (AIV09), and A/chicken/England/085598/2022 (AIV48). Infected epithelial cells were identified using immunohistochemistry in explants from both the teat and gland cistern for all three breeds following infection with AIV09 and AIV48, indicating that mammary tissue from each of the three tested cattle breeds cattle is permissive to H5N1 infection. Lectin staining showed expression of both 2,3-linked and 2,6-linked sialic acids in the mammary tissue of all donors showing that all three breeds have the potential to support infection with both avian-adapted and mammalian adapted IAVs. Together, these findings demonstrate that mammary glands from both beef and dairy cattle breeds are permissive to infection with avian-adapted and mammalian-adapted H5N1 viruses and highlight the potential for this tissue to act as a mixing vessel for IAV reassortment, underscoring the need to include cattle in ongoing H5N1 surveillance and risk-assessment frameworks. Impact StatementThe emergence of highly pathogenic avian influenza H5N1 in dairy cattle has expanded the recognised host range of influenza A viruses. Further, the ability of the virus to infect the mammary gland and transmit via milk revealed a novel interface for transmission to humans and animals. Although sustained circulation in US dairy herds has been reported, the susceptibility of mammary tissue from other breeds (including beef cattle) commonly used in different countries has been largely unexplored. Here, we show that avian-origin H5N1 viruses can infect tissues derived from the mammary gland of three common cattle breeds (Aberdeen Angus, Holstein Friesian, and Limousin). Virus was detected in epithelial cells from both dairy and beef breeds, indicating that H5N1 can infect multiple breeds. Receptor profiling showed abundant 2,3-linked and 2,6-linked sialic acids, consistent with a tissue environment that may support infection with both avian-adapted and mammalian-adapted viruses. These findings demonstrate that multiple cattle breeds are permissive to H5N1 infection and strengthens the evidence base for including cattle in H5N1 surveillance and risk-assessment frameworks.

Hepatitis D virus (HDV) is a satellite virus that utilizes hepatitis B virus (HBV) as a helper for infectious particle formation. HDV was originally identified as a novel antigen in liver biopsies of HBV patients, and later studies showed the "delta" antigen (DAg) to be the sole protein encoded by HDV. Until the discovery of HDV-like agents in birds and snakes in 2018, HDV was a unique example of animal satellite viruses. We identified Swiss snake colony virus 1 (SwSCV-1) in the brain of a Boa constrictor, and through comparison we found the genome organization of SwSCV-1 to resemble that of HDV. However, in addition to the DAg open reading frame (ORF), the genome of SwSCV-1 includes another >500 nt ORF, "ORF2". To study whether the putative ORF2-encoded protein plays a role in the SwSCV-1 life cycle, we established an infectious clone of the virus with a point mutation in the methionine initiation codon of ORF2. The mutation did not significantly affect initiation of replication, establishment of persistent infection, or infectious particle formation upon superinfection with a helper virus. Using additional methods, we gathered further evidence confirming that ORF2 is not actively translated in boa constrictor cells. We further showed that unlike HDV, SwSCV-1 expresses a single form of the DAg. Although the proteins encoded by SwSCV-1 and HDV only include one and two forms of the DAg, respectively, whether other kolmioviruses express additional forms of DAg or related proteins in some cell types or host species merits further research. IMPORTANCEApproximately 40 years after the discovery of hepatitis D virus (HDV), satellite viruses with similar genome organization were found in various animals, thereby giving rise to family Kolmioviridae. HDV encodes a single protein, the delta antigen (DAg), which comes in small and approximately 20 amino acids longer large form. The genome of some HDV species and many of the newly found kolmioviruses contains additional open reading frames (ORFs), potentially enabling protein expression. Here, we studied the viral proteins expressed during Swiss snake colony virus 1 (SwSCV-1) infection of boa constrictor cells. Our findings show that unlike HDV, SwSCV-1 encodes only a single form of DAg. In addition, our study suggests that, like in HDV, the additional ORF in SwSCV-1 genome does not give rise to a protein. Although we could not demonstrate expression of additional viral proteins during SwSCV-1 infection, it is important to study the proteome of other kolmioviruses.

Hantavirids, specifically the members within the genus Orthohantavirus, represent a significant global public health threat, with bat-associated lineages challenging traditional reservoir paradigms. To investigate the genetic diversity of hantavirids in Southeast Asia, we conducted an expanded surveillance program in Lao PDR from May 2023 to October 2025 in bat populations and wild animals from local wet markets. Using molecular screening and deep sequencing to characterize hantavirids from bat populations and wild animals from local wet markets, we identified 20 positive samples across four bat species, recovering coding-complete genomes for multiple novel variants. Phylogenetic analysis confirmed that these viruses form a monophyletic group within Mobatvirus, resolving into two major subclades. The first subclade clustered with Quezon and Robina viruses found in fruit-eating bats. The second subclade further split into two lineages corresponding to Thakrong and Xuan Son viruses, which are associated with trident and leaf-nosed bats, respectively. Despite the strong host specificity observed, the detection of these viruses in a wet market, a critical interface for human-wildlife contact, indicates a potential zoonotic risk. These findings significantly expand the known diversity of mobatviruses in Laos and highlight the urgent need for serological surveillance in at-risk human populations to assess the potential for spillover.

Since 2021, highly pathogenic avian influenza viruses (HPAIVs) belonging to H5N1 clade 2.3.4.4b have circulated widely in North American wild birds and repeatedly spilled over into mammals. In 2025, the first H5N1-associated deaths in humans were recorded in the Western hemisphere, raising questions about how the ongoing evolution of the virus in wild birds impacts spillover risk. Here, our analysis of 21,471 H5N1 genomes identified an evolutionary shift in mid-2024, driven by interhemispheric migration from Asia and reassortment with new antigens. The genotypes that dominated the early years of North Americas H5N1 epizootic traced their ancestry back to Europe, but Asia was the source of new "D1.1" genotype viruses that (a) spread faster, (b) have higher reassortment potential, (c) a broader host range, (d) repeatedly spill over to bovines, and (e) cause severe disease in humans, including non-farm workers.

The majority of emerging infectious diseases are zoonotic, having their origin in wildlife before spilling over into the human population. While small mammals are recognized as critical reservoirs for these viruses, their viral diversity remains largely uncharacterized across many African countries. We conducted molecular surveillance of synanthropic rodents and shrews in the Kibera informal settlement in Nairobi and the rural Taita Hills region of Kenya to detect and characterize potential zoonotic viruses. Tissue samples from 228 rodents and shrews were screened for six viral families using PCR assays. Rat hepatitis E virus (HEV) (Rocahepevirus ratti), a rodent-associated virus with potential for human spillover, was identified in Mus musculus and Rattus norvegicus from Kibera. NGS was conducted for the HEV positive samples, and we obtained two near-complete HEV genomes from Rattus norvegicus, which clustered within rodent-associated HEV genotypes in the phylogenetic analysis. The two sequences from the Rattus norvegicus cluster together, indicating a close genetic relationship. Paramyxoviruses belonging to the genera Jeilongvirus and Parahenipavirus were detected both from Taita and Kibera in nine different samples from Rattus norvegicus, Mus minutoides, Crocidura sp and Acomys ignitus. One paramyxovirus positive sample (Acomys ignitus) from Taita was selected for further sequencing with NGS, and a complete genome of a new jeilongvirus was assembled. Phylogenetic analysis of the detected viruses confirmed the close relation to previously known rodent-borne jeilongviruses but also revealed potentially novel jeilong- and parahenipavirus species. Our findings highlight the circulation of potentially zoonotic viruses in both urban and rural small mammals in Kenya. It emphasizes the necessity of continued genomic surveillance of zoonotic viruses to mitigate risks of their spillover into human populations. HighlightsO_LISurveillance reveals diverse rodent-borne viruses circulating in Kenya. C_LIO_LIRat-HEV was detected in Rattus norvegicus and Mus musculus from an urban low-income area. C_LIO_LIParamyxoviruses were detected across multiple rodent and shrew species, including novel Acomys ignitus jeilongvirus. C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=139 SRC="FIGDIR/small/719784v1_ufig1.gif" ALT="Figure 1"> View larger version (66K): org.highwire.dtl.DTLVardef@194e81eorg.highwire.dtl.DTLVardef@11342cdorg.highwire.dtl.DTLVardef@186ad97org.highwire.dtl.DTLVardef@eeb516_HPS_FORMAT_FIGEXP M_FIG C_FIG

Coat protein (CP) tertiary structures and their capsid organization of spherical viruses are highly conserved within each virus family. While AlphaFold successfully predicts the tertiary structures of individual CPs, their association to form proper quaternary assemblies cannot be easily accomplished. Here, we report a generalized methodology and associated web-based utility (https://foldavirus.org) that combines AlphaFold predictions of CPs with the knowledge on corresponding icosahedral architectures (e.g., T=1, 3, 4...) based on the known structures from the same virus family to generate associated capsids. The resulting assemblies are subjected to Amber energy minimization to relieve any steric clashes at the inter-subunit interfaces. Significantly, the capsid models are validated by calculating robust Mahalanobis distance using the residue annotations categorized as interface, core and surface amino acids with respect to those observed in the experimentally determined analogous structures. Given the amino acid sequence of CP(s), we successfully generated capsids up to T=9 icosahedral symmetry, including those of Picornaviruses that display pseudo-T=3 symmetry comprising different CPs. As the number of currently available CP sequences are 3-4 orders of magnitude larger than the experimentally determined 3D-structures, this approach bridges the huge gap that exists between the corresponding sequence and structure space.

Thogotoviruses are a group of tick-borne, six-segmented, negative-sense single-stranded RNA viruses. These viruses encode an RNA-dependent RNA polymerase that recognizes promoter sequences located at the genomic termini to initiate RNA synthesis. The 5' and 3' ends of the genome bind to the polymerase and function as a promoter. Outside the catalytic center, they base-pair with each other to form a double-stranded RNA structure. This structure is referred to as the distal duplex and plays an important role in RNA synthesis. In this study, we investigated how the RNA sequence of the distal duplex influences polymerase activity using minigenome systems of two thogotoviruses, Oz virus (OZV) and Dhori virus (DHOV). Each virus exhibits distinct activities among its six segments. In OZV, one determinant of these differences is the base pair at positions 5'12 and 3'11 within the distal duplex, where promoter activity varies depending on whether the base pair is G:C or A:U. In contrast, the DHOV polymerase is not affected by this difference. These results indicate that, even within the genus Thogotovirus, viruses differ in whether they possess a mechanism that modulates promoter activity based on subtle sequence differences within the distal duplex. Furthermore, phylogenetic analysis and comparison of promoter sequences suggest that thogotoviruses can be divided into groups that do or do not regulate intersegment promoter activity via the base pair at positions 5'12 and 3'11. HighlightsO_LIMinigenome systems of Oz virus and Dhori virus reveal segment-specific differences in promoter activity C_LIO_LIThe distal duplex sequence modulates RNA synthesis in a virus-dependent manner C_LIO_LIThe base pair at positions 5'12/3'11 determines promoter activity in Oz virus but not in Dhori virus C_LIO_LIThogotoviruses can be divided into groups that do or do not regulate promoter activity via distal duplex sequence variation at positions 5'12/3'11 C_LI